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There exist important stochastic physical processes involving infinite mean waiting times. The 
mean divergence has dramatic consequences on the process dynamics. Fractal time random walks, 
a diffusion process, and subrecoil laser cooling, a concentration process, are two such processes 
that look qualitatively dissimilar. Yet, a unifying treatment of these two processes, which is the 
topic of this pedagogic paper, can be developed by combining renewal theory with the generalized 
central limit theorem. This approach enables to derive without technical difficulties the key physical 
properties and it emphasizes the role of the behaviour of sums with infinite means. 
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Introduction 

The fractal time random walk 0,0 has been developed 
in the 1970's to explain anomalous transport of charge 
carriers in disordered solids. It describes a process in 
which particles jump from trap to trap as a result of 
thermal activation with a very broad (infinite mean) dis- 
tribution of trapping times. It results in an unusual time- 
dependence of the position distribution which broadens 
while the peak remains at the origin. The method of 
choice to study the fractal time random walk is the con- 
tinuous time random walk technique. 

Subrecoil laser cooling 0, Q has been developed in the 
1990's as a way to reduce the thermal momentum spread 
of atomic gases thanks to momentum exchanges between 
atoms and laser photons. It is a process in which, as 
a result of photon scattering, atoms jump from a mo- 
mentum to another one with a very broad distribution of 
waiting times between two scattering events. It results 
in an unusual time dependence of the momentum distri- 
bution which narrows without fundamental limits hence 
giving access to temperatures in the nanokelvin range. 
The method of choice to study subrecoil laser cooling is 
renewal theory 

Fractal time random walks and subrecoil cooling seem 
at first sight very dissimilar. The first mechanism gener- 
ates a broader and broader distribution, while the second 
generates a narrower and narrower distribution. Nev- 
ertheless, inspection of the theories of both phenomena 
reveals strong similarities: the continuous time random 
walk and the renewal theory are two closely related ways 
to tackle related stochastic processes. Physically, the two 
mechanisms share a common core, a jump process with 
a broad distribution of waiting times. 

The aim of this pedagogic paper is to bridge the gap 
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between fractal time random walk and subrecoil laser 
cooling. We show that the essential results of the two 
theories can be obtained nearly without calculation by 
combining the simple probabilistic reasoning underlying 
renewal theory and the generalized central limit theorem 
applying to broad distributions. This provides more di- 
rect derivations than in earlier approaches, at least for 
the basic cases considered here. 

In the first part, we describe the microscopic stochastic 
mechanisms at work in the fractal time random walk and 
in subrecoil cooling and relate them to renewal processes. 
In the second part, we explain elementary properties of 
renewal theory and derive asymptotic results using the 
generalized central limit theorem and Levy stable distri- 
butions. In the third part, we draw the consequences for 
the fractal time random walk and subrecoil cooling. The 
fourth part contains bibliographical notes. 



I. FRACTAL TIME RANDOM WALK AND 
SUBRECOIL COOLING: MICROSCOPIC 
MECHANISMS 

A. Fractal time random walk 

The notion of fractal time random walk emerged from 
the observation of unusual time dependences in photo- 
conductivity transient currents flowing through amor- 
phous samples. It can be schematized in the following 
way. 

Consider first a one dimensional situation called the 
Arrhenius cascade Q in which the charge carriers are 
placed in a random potential with many local wells and 
barriers and can jump from one well to another one 
thanks to thermal activation (Fig. QJl). The Arrhenius 
cascade potential presents two features: a global tilt rep- 
resenting the effect of the electric field on the carriers and 
local random oscillations creating metastable traps sep- 
arated by barriers representing the disorder created by 
the amorphous material. Thus the potential seen by the 
carriers is a kind of random washboard with a discrete 
number of metastables states. 

The mean lifetime of state i, i.e., the mean waiting 
time before the occurrence of a thermal jump, is given 
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FIG. 1: Biased random walks in a disordered system (a) Ar- 
rhenius cascade. Each carrier is put in random potential with 
a global tilt. It undergoes jumps over barriers of random 
heights Ei from a metastable well to the next one on the right 
thanks to thermal activation, (b) Photoconductivity setup. At 
time t = 0, a light pulse creates carriers in the immediate 
vicinity of the left electrode of an amorphous sample. The 
carriers of one sign then move through the sample thanks to 
an applied electric field. (The carriers of the opposite sign are 
immediately absorbed by the left electrode.) 



by the Arrhcnius law: 



= r e^ kT 



(1) 



where To is a time scale, is the height of the energy 
barrier separating state i from state i + k is the Boltz- 
mann constant and T is the temperature. The potential 
global tilt is assumed to be large enough to neglect back- 
ward jumps from i to i— 1. The random walk we consider 
is thus completely biased. The time spent between the 
metastable states is neglected. For a given barrier height 
Ei, the lifetime distribution xj)(j\Ei) is exponential with 
mean f, : 



i>{r\Ei) = - e 

T, 



— t/t 



(2) 



barrier height Ei (eq. J5J) with the barrier height statis- 
tics conveniently described by an exponential distribu- 
tion P(Ei), 



P(E; 



i) = 7^ e~ E ' /E ° for E t > 0, 



E, 



(3) 



where Eq is an energy scale related to the sample disor- 
der. 

The waiting time distribution ip( T ) is then 

poo a 

1>(t) = J dE t PiE^irlEi) = 07(1 + a, r/r ) -g- 

T (4) 

where j(a',x) = J Q e~ u u a _1 du is the incomplete 
gamma function and 
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At long times, i(j(t) tends to a power law, hence the term 
"fractal time random walk" : 



^(r)ca 2 r(a)^ 



(6) 



with r(a) = / °° ii -^-" du. 

If a < 1, states i have an infinite mean lifetime 
(r) = J Q Tip(r) dr. However, they are unstable since 
they all ultimately decay to the next state + Usu- 
ally, unstable states have a well defined and finite mean 
lifetime. Here, the somewhat paradoxical presence of un- 
stable states with infinite mean lifetimes is at the origin of 
the striking properties of the fractal time random walk. 



B. Subrecoil laser cooling 



In the photoconductivity experiments (Fig.^>), a light 
pulse creates at time t = carriers localized near the 
surface of the sample. The carriers then move through 
the disordered sample thanks to an applied electric field. 
Thus, this situation can be modelled by a large number 
of Arrhenius cascade in parallel, each electron path being 
associated to one cascade. 

One may (wrongly) expect that the transient current 
flowing through the sample is quasi-constant at the be- 
ginning, while the bunch of carriers propagates through 
the sample, before decreasing rapidly to zero when the 
carriers leave the sample after reaching the end electrode. 
But what is observed is quite different. The current de- 
creases as a power law ~ l/t l ~ a while the carriers are 
still in the sample, then as ~ l/t 1+a when some carriers 
start leaving the sample. For simplicity, we assume here 
that the sample is semi-infinite so that the carriers never 
leave the sample. 

The explanation of this anomalous behaviour will be 
shown to be related to the distribution of lifetimes r,. 
The randomness of the Tj 's results from the combination 
of the exponential statistics of jump times for a given 



Laser cooling of atomic gases consists in reducing the 
momentum spread of atoms thanks to momentum ex- 
changes between atoms and photons. Subrecoil laser cool- 
ing consists in reducing the momentum spread to less 
than a single photon momentum, denoted hk. This para- 
doxical goal is achieved by introducing a momentum de- 
pendence in the photon scattering rate (see Fig. so 
that it decreases strongly or even vanishes in the vicinity 
of p = 0, where p denotes the atomic momentum, taken 
in one dimension for simplicity. 

The mechanism of subrecoil cooling is explained in 
Fig. [3 Any time a photon is absorbed and spontaneously 
reemitted by an atom, the atomic momentum undergoes 
a momentum kick on the order of hk, which has a ran- 
dom component because spontaneous emission occurs in 
a random direction. Thus, the repetition of absorption- 
spontaneous emission cycles generates for the atom a mo- 
mentum random walk (see Fig. |2t>) , with momentum de- 
pendent waiting times r between two kicks. When an 
atom reaches by chance the vicinity of p = 0, it tends to 
stay there a long time. This enables to accumulate atoms 
at small momenta, i.e., to cool. 
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FIG. 2: Subrecoil laser cooling (a) The mean sojourn time at 
momentum p f(p) becomes very large for small atomic mo- 
menta, (b) Photon scattering creates a momentum random 
walk with an accumulation in the vicinity of p = due to the 
momentum dependence of the mean sojourn time f(p). 



For a quantitative treatment, we introduce f (p), the 
mean sojourn time at momentum p (also the mean wait- 
ing time between two spontaneous photons for an atom 
at momentum p). For a given p, the sojourn time at mo- 
mentum p, i.e., the distribution ip(r\p) of sojourn times 
at momentum p is 



iP( t \p) 
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We need to characterize the distribution of "landing" 
momenta 7r(p) after a spontaneous emission. Under 
favourable but often realistic assumptions, atoms spend 
most of the time around the origin in the interval 
[— ptrap * +Ptrap] > because they diffuse fast outside this in- 
terval and thus come back to it rapidly after leaving it. If 
Ptrap < hk, then after a spontaneous emission, the distri- 
bution of atomic momenta can be considered as uniform: 



7r(p) 



(8) 



The distribution of sojourn times after a spontaneous 
emission is thus 



Ptrap 



dp ir(jp)ip(r\p). 



(9) 
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We consider the physically relevant case of power law 
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where (3 > 0, and r and po are time and momentum 
scales, respectively. Then, one finds, just as in the fractal 



time random walk, a waiting time distribution with a 
power law tail: 
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If P < 1, the mean waiting time is finite and simple 
integration gives 
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If j3 > 1, on the contrary, the mean waiting time is in- 
finite. This divergence of the mean has dramatic (and 
positive in terms of cooling) consequences (see flIII Bp . 



C. Connection with renewal theory 

Renewal processes are stochastic process in which a 
system undergoes a sequence of events (denoted by • in 
Fig.yJ) separated by independent random "waiting times" 
ti, T2, ... The term "renewal process" comes from en- 
gineering. Assume that, at time t = 0, one installs a 
machine in a factory. When, after being operated for a 
random lifetime n, the machine breaks down, it has to 
be replaced by a new one, which will work till it breaks 
down at t± +T2 and has to be replaced ... If, instead of a 
single machine, one has installed a large number of iden- 
tical machines, then, to decide how many replacement 
machines must be stored at a given time, one needs to 
know the replacement rate, which we call hereafter the 
renewal density. 

< X > 
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FIG. 3: Renewal processes. The system undergoes a sequence 
of events (jumps from trap to trap, momentum kicks ...) at 
random times separated by waiting times n, T2, ... 

To understand the statistical properties of renewal pro- 
cesses, various quantities are introduced. The most de- 
tailed information is provided by the distribution of the 
number of renewals, ft(r), i.e., the probability distribu- 
tion for the system to undergo r events in time t. One also 
introduces derived quantities, the mean number of re- 
newals at time t, (r) t , and the mean renewal rate at time 
t, denoted R(t) and called the renewal density. Math- 
ematical expressions for these three quantities will be 
given in £111 Al Here we show the role they play in fractal 
time random walk and in subrecoil laser cooling. 
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In the biased (fractal time or non fractal time) random 
walk, the discretized positions n at time t correspond di- 
rectly to the number of jumps performed between time 
and t. Hence the position distribution p(n, t) in the frac- 
tal time random walk is the renewal number distribution: 



p(n,t) = f t (n). 



(14) 



The mean position of the carriers is (r) t - The current 
i(t) measured in photoconductivity experiments before 
the carriers get out of the sample is proportional to the 
mean carrier velocity d(r) t /dt, which is the renewal den- 
sity (see dim . Thus, one has 



i(t) (X R(t). 



(15) 



In subrecoil cooling, the momentum distribution p(p, t) 
can be written in the following form: 



p{p,t) = tt(p) / dti R(ti)P a (t - t t \p), 



(16) 



where ti is the time of the last jump occurring between 
and t, R(ti)dti is the probability that a jump occurs 
during the interval [U,ti +d£/), 7r(p) is the uniform prob- 
ability distribution for a jumping atom to land at mo- 
mentum p and P s (t — ti \p) is the survival probability for 
an atom landing at momentum p at time U to stay there 
till at least time t. Using eq. Q, one has trivially 



P s (t-ti\p) 



dr ip{ T \p) = e 



- e -(t-t,)/f(p) 



(17) 



The non trivial physical information is contained in the 
renewal density R(t). 

The height h(t) of the momentum distribution peak, 



h(t)=p(p = 0,t), 



(18) 



is proportional to (r) t , the mean number of jumps be- 
tween and t. Indeed, for any jump, there is a proba- 
bility 7r(p)2dp = dp/ptrap to fall in the vicinity [—dp, dp] 
of the origin and to stay there indefinitely since states 
in [—dp, dp] have arbitrarily long lifetimes in the limit 
dp — * (f(p) — * oo, see eq. {TO}). Thus the height 

writes 



p-*0 



h(t) 



2p 



'trap 



(19) 



II. RENEWAL THEORY AND LEVY STABLE 
LAWS 

A. General formulae 

The number of renewals r t in a time t is defined as 
the number of jumps having occurred before time t. It 
satisfies 



S rt < t < Sy t + 1 



(20) 



where S rt = Y^iLi T i 1S the sum of the first r t waiting 
times. The relationship between the renewal number dis- 
tribution f t (r) and the waiting time distribution ^(r) can 
be obtained from the following simple reasoning. 

Note first that the distribution, denoted ip r *(S r ), of 
the sum S r of r independent identically distributed wait- 
ing times is the r th convolution product of ip(r) with 
itself. Moreover, from the definition of the number r t of 
renewals, one has obviously 



Pr(r t < r) = Pr(S r >«) = [!- * r *(t)] 



(21) 



where 1 i> r *(S r ) — J Q r ip r *(u) du denotes the distribution 
function of S r (in spite of its notation, ty r *(S r ) is not the 
r th convolution product of the waiting time distribution 
function ^(t) = J Q i[>(u) du). The probability distribu- 
tion /t(r) of the number of renewals at time t is thus 
finally 

f t (r t = r) = Pr(r t < r+l)-Pr(r t < r) = * r *(t)-* (r+1) *(t). 

(22) 

This expression relates the distribution of a discrete ran- 
dom variable, r, to the distribution fonctions of contin- 
uous random variables t. Important quantities derived 
from the renewal number distribution /t(r) are the mean 
number of renewals at time t: 



( r )t =22rf t {r) 

r=0 

and the renewal density, i.e. 
newals per unit time|l9|: 



(23) 



the mean number of re- 



R(t) 



d(r) 
df 



(24) 



Usual theoretical treatments of renewal problems are 
based on Laplace transforms of the waiting time distribu- 
tion that are indeed well suited to handle the convolution 
products ip r *(t). Here, we prefer a different approach 
based only on the generalized central limit theorem. This 
approach requires nearly no calculation. Moreover, it 
stresses the role of the behaviour of sums of many ran- 
dom variables. 

Indeed, most useful distributions tend, under repeated 
convolution, to a Levy stable law given by the generalized 
central limit theorem. Thus we can obtain from eq. ill' 21) 
analytical expressions for f t (r) and for related quantities 
in the limit of large r and hence large t. Depending on 
the finitcness of the first two moments of the waiting time 
distribution ip(r), three cases can be distinguished. For 
brevity, we treat only the two most striking cases: first 
and second moment both finite in 311 Bl first and second 
moment both infinite in 311 CI 



B. Case of waiting time distributions with finite 
first and second moment 



When both the mean p = (r) and the variance a 2 of 
the waiting time are finite, the (usual) central limit the- 
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orem applies and ip r * (t) tends to a Gaussian distribution 



■ exp 



2ol 



(25) 



with mean \i r = ru and variance of = rcr 2 . Thus, after 
changing variables in eq. (|22|l . one has 



M r ) 



('-Mr+l)K+l 



2 / 2 du. 



(26) 



Fix (t — n r ) I 'ay and take t — > oo. One expects intuitively 
that the number of renewals is approximately given by 
t/fx ^> 1 at large t (this is validated a posteriori by 
eq. pgfr). Thus (t - Mr)/or - (< - fir) / (a^/t/fi) and 

(t-/x r+ i)/(T r+ i ~ (t - fir) / (a y/t/Ji) - This leads 

for the number of renewals to a Gaussian distribution 
with mean t/ a and variance a 2 t/ fi 3 : 



Mr) 



■ V /27Tg 2 t/^ 3 



exp 



(r-t/M) 2 

' 2(T 2 t//i3 



(27) 



Although well known in renewal theory, this result is non 
trivial: the distribution of the sums of r (— > oo) terms 
and of the number of renewals at large times are found to 
have the same (Gaussian) shape. This is grossly violated 
with infinite mean waiting times (see flll C() . 

As direct consequences of eq. I|27|l. the mean number 
of renewals (r)t tends to t/fi at large times: 



<r>i 



t 



(28) 



and, using eq. (|24p. the renewal density tends to the re- 
ciprocal of the mean waiting time: 



R(t) 



1 

/' 



(29) 



in agreement with intuition. 

Physical comments and an example of a renewal pro- 
cess with finite mean waiting time will be presented in 

mm 



C. Case of waiting time distributions with infinite 
first and second moment 



We consider now the case of waiting time distributions 
with infinite first two moments, focusing on the canonical 
example of distributions with (Pareto) power law tails of 
index a, 



V(r) 



T-l + Q : 



(30) 



with < a < 1 (slightly more general cases can be han- 
dled using the theory of regular variations). The distri- 
butions ip r * of the sums S r no longer tends to Gaussians 



but, according to the generalized central limit theorem, 
to Levy stable laws: 



V*(t) 



„l/a 



Ln 



„l/a 



(31) 



where L _b(u) is a one-sided (it > 0) Levy stable law 
given by its Laplace transform, 



I L Qj s(u) du = e 



-Bs 



(32) 



and B is a scale parameter given by 

B = r(l-a)T Q . (33) 
After changing variables in eq. (|22|l , one obtains 

Mr) = / L a , B (u) du, (34) 

Jt/(r+iy/» 

which leads to the following asymptotic expression 20] : 
, m .i+i/ Q L*,b (pys-) (35) 



Mr) 



In sharp contrast with the finite moments case, the re- 
newal number distribution f t (r) differs strongly from the 
distribution of the sum (eq. (|31|l ) even though both dis- 
tributions involve Levy stable laws. As will be shown in 
fllll Al f t (r) has a slow decay at small r and a fast decay 
at large r unlike -0"(i) which decays as l/t 1+a . 

The mean number of renewals (r) t is finite and can be 
related to a negative moment of a Levy stable law: 



/>oo 

a / u~ a L a , B {u) du. (36) 



rf t (r) dr = t 



Hence, using f °° u a L a ,B(u) du — l/(BaT(a)) and 
cq. I|33|l . one has finally 



sin(7ra) 



(37) 



Finally, using eq. (|24|) , the renewal density tends asymp- 
totically to an ever decreasing power law: 



R{t) 



sin(7ra) 



1 



— a/1- a ' 



(38) 



Physical comments and an example of a renewal pro- 
cess with infinite mean waiting time will be presented in 

mm 



III. APPLICATION TO FRACTAL TIME 
RANDOM WALK AND SUBRECOIL COOLING 

We have seen in ill CI that the most important quan- 
tities appearing in biased random walks and in subre- 
coil cooling are directly related to renewal theory: re- 
newal number distribution ft(r) (position distribution 
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in biased random walks), mean number of renewals (r)t 
(peak height in subrecoil cooling) and renewal density 
R(t) (current in the biased random walk, momentum dis- 
tribution in subrecoil cooling). 

In this section, we thus apply the results on renewal 
processes obtained in SjH]to biased random walks and sub- 
recoil cooling. We emphasize the physics consequences 
of the divergence of the mean waiting time in renewal 
processes by opposing, for each application, one exam- 
ple with finite mean waiting time and one example with 
infinite mean waiting time. These examples reveal the 
generic features of the cases with finite or infinite mean 
waiting times. 

A. Biased random walks 
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As an example of waiting time distribution with finite 
mean and standard deviation generating a (non fractal 
time) random walk, we consider an exponential distribu- 
tion 



V>(r) = - e~ T ^ for r > 
A* 



(39) 



with mean and standard deviation both equal to p. In 
this case, the convolutions of ip(r) have the simple ex- 
plicit forms of the Gamma (or Erlang) distributions 



(n- 1)!// T| 



-t/h 



(40) 



Thus, applying eq. (J22J, the renewal density f t (r) and 
hence the position distribution p(n, t) are exactly known: 



p{n,t) 



-t/fi 



(41) 



(One recognizes the Poisson distribution of mean (r)t = 
t/n, as expected: in this case, the renewal process is a 
Poisson process.) As shown in Fig.QJ this exact position 
distribution rapidly tends to the Gaussian distribution 
given by eq. (|2"7|l : 



p{n,t) 



v / 2nt/Jl 



exp 



It III 



(42) 



We thus recover the intuitive picture of normal trans- 
port: a position distribution propagating at a constant 
speed 1/p and spreading as \ft. The current i(t) result- 
ing from such a distribution in photoconductivity exper- 
iments is related to the renewal density (eq. I|15[)) which 
leads, using eq. (|2UJ), to 



i{t) (x 



(43) 



As an example of waiting time distribution with in- 
finite mean and standard deviation generating a fractal 



FIG. 4: Time evolution of the position distribution p(n, i) for 
a biased random walk with an exponential waiting time distri- 
bution [finite (t)). The distribution propagates at constant 
speed 1/u and spreads as \ft. The exact result (eq. H4H - ) and 
the asymptotic one (eq. (1421 s ) agree even at short times. 



time random walk, we consider a Pareto distribution |2] 
of index a = 1/2, 



V(t) 



1/2 
2r 3 / 2 



for r > tq 



(44) 



with To > 0. As there is no simple analytic form for 
ip n *(t) and thus for the position distribution p(n,t) = 
# n *(i) - ( eq . (j22J)), we obtain p(n,t) through 

numerical simulation (see Fig. [SJ. On the other hand, 
using eq. I|35|) and the known form for the asymmetric 
Levy stable law of index 1/2 called the Smirnov law (or 
using eq. H32[) h the asymptotic distribution is found to 
be a half-Gaussian 



p{n,t) 



[TQ ( 7TT 

— exp r 



(45) 



The fact that a half-Gaussian is obtained should not give 
the impression that fractal time transport is similar to 
normal transport described by full Gaussians (eq. H42|) h 
First, the half-Gaussian is specific to Pareto waiting time 
distributions with a — 1/2 (see below for other a's). Sec- 
ond, the properties of the half-Gaussian in fractal time 
transport are completely different from those of the full 
Gaussian of the normal transport. Indeed, instead of 
propagating, the distribution peak remains at the origin 
n = at all times. Only the tails spread to the right, 
more and more slowly as times goes by. This is due to 
the fact that the carriers statistically tend to be trapped 
into deeper and deeper traps at long times, which slows 
down their motion. The resulting current is thus a de- 
creasing function of time, at all times. Using eq. I|15|) and 
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eq. one obtains: 



i(t) 



(46) 



up to power law corrections. As 1/(1 — a) > 1, the trans- 
port front always decreases faster than exponentially and 
presents a well defined characteristic position. Conse- 
quently, the mean carrier position yielding the current, 
is well defined, unlike the mean waiting time which is 
infinite. Using eq. Ijl5|) and eq. (|38|1 . one finds that the 
current decays as a power law: 
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FIG. 5: Time evolution of the position distribution p(n, i) for 
a biased fractal time random walk constructed from a Pareto 
distribution of index a = 1/2 {infinite (t)). The distribution 
spreads slowly towards positive values but its maximum al- 
ways remains at the origin. Note the smaller position scale 
and the longer time scales compared to Fig. |1] which em- 
phasizes the slowness of transport in the fractal time random 
walk. The exact (simulated) and asymptotic results (eq. 1)45^ ) 
are in good agreement. (They have seemingly different norms 
at short times because the exact result is a discrete distribu- 
tion while the asymptotic one is a continuous distribution.) 

It is worth examining the general case of waiting time 
distributions with power law tails like eq. Ij30|l and in- 
finite means (a < 1). From the following expansion of 
asymmetric Levy stable laws with a < 1, 



aB 



(47) 



*oo r(l - ajxi+'x ' 
one finds the behaviour close to the origin (for t — + oo): 

P M n% ^y. (48) 

Hence, the position distribution is flat (n independent) at 
small n, as already observed for the special case a = 1/2. 
Moreover, using 



L aB (x) ~ Ac 2 ^-") exp 



C 



x a/(l-u) 

where A and C are constants, one obtains 

C r l/(l-a) 



p(n, t) cx exp 



(49) 



(50) 



lit) OC -: 

v ' 4-1 — a 



(51) 



as already observed for the special case a — 1/2 and 
in agreement with photoconductivity transient experi- 
ments. 



B. Subrecoil cooling 

Consider first the case of waiting time distributions 
with a — 1/P > 2 (eq. ensuring a finite mean wait- 

ing time (r) (and a finite (r 2 )). According to eq. f2"5|l 
and eq. (I3)l. one has 



R(t) 



^ 1 - §_ ( Ptrap 



(52) 



and thus, applying eq. l(TfJ|) with eq. |JS} and eq. (|T7|) . one 

finds the momentum distribution 



1-/3 



1 



Ptr 



IpI 







2.Ptrap 

This distribution has stationary tails: 



1/73 



(53) 



(54) 



but a non stationary peak that increases linearly in time: 



pip, t) 



1-/3 



— for \p\ < po 

Tn \ t 



1//3 



(55) 



This peak is also obtained directly using relation l|19|) 
between the height h(t) and the mean number of renewals 
of eq. (J2EI- 

Finite mean waiting times are not very favourable for 
the cooling since only a vanishingly small fraction of 
atoms goes on accumulating at smaller and smaller ve- 
locities (see Fig. [HJ - 

Consider now the case of waiting time distributions 
with infinite mean waiting times (a = 1/0 < 1) ensuring. 
According to eq. i|38[l . one has 



R(t) 



sin(7ra)p t ra 



1 



t^oo TraT(a)po T^t 1 a 



(56) 



(one must write ip(t) of eq. in the form of eq. 1301 1. 
thus replacing by aT(a){pa/p tTa , v )Tg in eq. ggj). 
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FIG. 6: Time evolution of the momentum distribution p(p,t) 
for subrecoil cooling with a finite mean waiting time. Param- 
eters: P = 0.25, po = 1. The tails reach a stationary state. 
Only a vanishingly narrow part of the peak goes on increasing 
at long times (note the logarithmic p scale) . 



Thus, applying eq. (|16[) with eq. JSJ and eq. I|17|) . the 

momentum distribution writes 



P(p, t) 



where 



sin(-7ra) 
2ira 2 T(a)po 



Q(q) = a du 
Jo 



tp 



13 



t a-i e -(i-«)g 



(57) 



(58) 



is a confluent hypergeometric function. Thus p(p, t) 
presents a scaling form and evolves at all times scales, 
with no stationary state. The tails behave as 



p(p, t) 



(59) 

and the peak, also obtained directly from eq. (|19|l and 
the mean number of renewals (|37|) . behaves as 



P(P, t) 



sin(7ra) 
2Tta 2 r(a)pQ 



for 



H<»(?) 



1//3 



( 6 °) 

Infinite mean waiting times are favourable for the cooling 
(Fig. |7|): all atoms accumulate in a narrower and nar- 
rower peak in the vicinity of p = 0. The cooling goes on 
without fundamental limits at long times. The absence 
of limits is related to the significant weight (c/. (r) = oo) 
of p states with arbitrarily long lifetimes. 



IV. BIBLIOGRAPHICAL NOTES 

Fractal time The theory of fractal time random walks 
was developed in particular cases in [?J using continu- 



FIG. 7: Time evolution of the momentum distribution p(p, t) 
for subrecoil cooling with an infinite mean waiting time. Pa- 
rameters: /3 — 2, po — 1. The full momentum distribution, 
tails and peak, never reaches a stationary state. All atoms 
go on accumulating without limit in a narrower and narrower 
peak as time increases. 



ous time random walks and, in the general case, in 0, 
adding Tauberian techniques to handle Laplace trans- 
forms. Fractal times were used in Q to model turbu- 
lent diffusion. Several other applications involving fractal 
times were presented in and 0- 

Subrecoil cooling The statistical approach to subrecoil 
cooling was developed in three steps: 01 > E3 an d IH- 
The connection of this approach with renewal theory was 
stressed in 

Renewal processes The basic theory of renewal pro- 
cesses is presented in a small book [l^. This book does 
not include infinite mean waiting times but some general 
. ^expressions for finite mean waiting times are still valid in 
' this case too (see flll A|) . Ref. [Lj, chapter XI, presents 
renewal theory with a more theoretical viewpoint and in- 
cludes some results on infinite mean waiting times which 
were discovered at the beginning of the sixties, in spite of 
the fact that no application seemed to be known at the 
time. Renewal processes in a physics context including 
infinite mean waiting times have been studied systemat- 
ically in |14| using Laplace transform techniques. 

The mean number of renewals given by eq. <|37|l for 
infinite mean waiting times agrees with eq. (8) in 
obtained through Tauberian theorems for the fractal time 
random walk (there is to be a misprint in this reference: 
the expression for ip{r), page 424, line 3, must be replaced 
by iP(t) ~ a/ [t 1+a T(l + a)A(t)Vl It also agrees with 
the renewal theory developed in [lj] (eq. (3.6)). 

The renewal number distribution expressed with Levy 
stable laws given by eq. (|3*5|l agrees with eq. (5.6) of |13j . 
which is however more complicated. 

Generalized central limit theorem The generalized 
central limit theorem and results on Levy stable laws 
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used in £111 CI can be found, e.g. in Appendix B of [l^ 
or in §4.2 of Details on Levy stable laws can also be 
found in §1.2 of [3. 

Conclusion 

This paper has underlined the common statistical core 
at work in two seemingly opposite problems: a diffu- 
sion mechanism (fractal time random walk) and a cooling 
mechanism (subrecoil laser cooling). This core is made 
of waiting time distributions with infinite means. 

Usual theoretical techniques for these problems are 
based on Tauberian theorems which, through the Laplace 
transform, imply some loss of physical intuition. Here, we 



have presented a renewal theory approach which, thanks 
to the generalized central limit theorem, provides a short- 
cut to obtain physically relevant quantities. It empha- 
sizes the key contribution of the unusual behaviour of 
sums of random variables with infinite means. 

Understanding these stochastic processes with infinite 
means is not a purely academic game. It has already led 
to significant improvements of laser cooling strategies [l7| 
and other improvements are under way |18| . 
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